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Abstract 

Let the function u satisfy Dirichlet boundary conditions on a bounded 
domain SI. What happens to the critical set of the Ambrosetti-Prodi operator 
F{u) = —Au — f{u) if the nonlinearity is only a Lipschitz map? It turns out 
that many properties which hold in the smooth case are preserved, despite 
of the fact that F is not even differentiable at some points. In particular, a 
global Lyapunov-Schmidt decomposition of great convenience for numerical 
solution of F(u) = <? is still available. 

Keywords: Semilinear elliptic equations, Ambrosetti-Prodi theorem. 

MSC-class: 35B32, 35J9I, 65N30. 

1 Introduction 

A familiar set of hypotheses for the celebrated Ambrosetti-Prodi theorem is the 
following. Let P C K" be an open, bounded, connected domain with smooth 
boundary dD, and denote by 0 < Ai < A 2 < ... the eigenvalues of the free Dirichlet 
Laplacian —A on D. Let / : R —>■ K be a smooth, strictly convex function, with 
asymptotically linear derivative so that 

Ran /' = (a, &), a < Ai < & < A 2 . 

Under such hypotheses, the theorem states that the equation 

F{u) =-Au-f{u) = g, u|an = 0 (1) 

for, say, g S has (exactly) zero, one or two solutions in 

1.1 A first approach — locating C and F{C) 

In the original arguments m. [H]), a fundamental role is played by the critical 
set C of A : X = Y = Here is the subspace of 

functions of satisfying Dirichlet boundary conditions. The proof follows a 

few steps, which follow from the inverse function theory and the characterization of 
fold points. 

1. The critical set C C A is a hypersurface; every critical point is a fold. 
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2. F is proper, its restriction to C is injective and F ^{F{C)) = C. 

3. The spaces X — C and Y — F{C) have two components. Each component of 
X — C is taken injectively to the same component of E — F{C). 

1.2 A global Lyapunov-Schmidt decomposition 

Berger and Podolak came up with a different approach [2] , which is easier to phrase 
ioT F : X ^ Y between Sobolev spaces, X = and Y = ~ F{q{Q). 

Their main result is the construction of a global Lyapunov-Schmidt decomposition 
for F. Let ipi be the (positively normalized) eigenfunction associated to the ground 
state Ai, set Vx =Vy = (yi) and consider the orthogonal decompositions 

X = Vx®Wx and Y = Vy®Wy 

into vertical and horizontal subspaces. With a different vocabulary, their proof 
essentially goes through the verification of the following properties, by making use 
of spectral estimates on the Jacobians DF{u) : X —^ T. 

1. Horizontal affine subspaces of X are taken by F to sheets. 

2. The inverse under F of vertical affine subspaces of Y are fibers. 

3. Sheets are essentially flat, fibers are essentially steep. 

An affine horizontal subspace of X is a set of the form x + Wx , for a fixed x G X-, 
an affine vertical subspace of Y is of the form y + Vy, tor y G Y. Sheets are graphs 
of smooth functions from Wy to Vy and fibers are graphs of smooth functions from 
Vx to Wx. The third property states that the inclination of the tangent spaces to 
sheets, with respect to horizontal subspaces, and to fibers, with respect to vertical 
subspaces, is uniformly bounded from above. 

In the words of [4], F is a flat map. In the Ambrosetti-Prodi case, vertical 
spaces and fibers are one dimensional. More generally, the dimension equals the 
number k of eigenvalues of — in the set /'(K) (we suppose non-resonance, i.e., 
the asymptotic values of Ran/' are not eigenvalues) and similar results still hold. 

1.3 Sheets and fibers allow for weaker hypothesis 

There is an interesting bonus obtained from considering this global Lyapunov- 
Schmidt decomposition. For many nonlinearities /, the related nonlinear map 

F :X ^Y, Fiu) = -Au - f{u) 

is not proper, and part of the technology related to degree theory simply breaks 
down. The Ambrosetti-Prodi hypothesis yield properness of F, but it is not really 
essential to most of what we want to do. The first section of the paper is dedicated 
to explicit examples of Lipschitz nonlinearities for which properness does not occur 
— there are functions g GY for which F~^{g) contains a full half-line of functions 
in X. From such examples, we may obtain smooth nonlinearities with the same 
property, but we give no details. 

The results in the paper indicate that in most directions, properness holds. 
What we mean by this is something similar to the fact that even when a map 
F does not have an invertible Jacobian DF{u) at a point u, it may still have a 
subspace L{u) on which the restriction of DFiu) acts injectively — this is how 
bifurcation equations come up: they concentrate in a possibly small subspace S{u) 
transversal to L[u) the difficulties which are not resolved by the linearization at 
u. What we shall see is that the domain X of F splits into special (nonlinear) 
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surfaces of finite dimension k, on which properness may break down, which have for 
tangent spaces the subspaces S{u) when u is critical, but on transversal directions 
to these surfaces, the horizontal affine subspaces, properness is always present. As 
we shall see, counting or computing solutions of the differential equation F[u) = g 
boils down to a finite dimensional problem, simplifying both (abstract) analysis and 
numerics. 

1.4 A related numerical algorithm to solve the PDE F{u) = g 

Smiley and Chun ([l2],[T3]) showed that an analogous global Lyapunov-Schmidt 
decomposition exists for appropriate non-autonomous nonlinearities f{x,u{x)) and 
emphasized its relevance for numerical analysis: one might solve F{u) = g hy 
restricting F to the fiber ag containing F~^{g). In [4], this project is accomplished: 
given a right hand side g, one first obtains numerically a point in ag, which in 
the Ambrosetti-Prodi case is a curve, and then proceeds to search for solutions 
by moving along it. The algorithm is sufficiently robust to handle more flexible 
nonlinearities: it does not require convexity of / or properness of F, and the range 
of /' may include other (finite) sets of eigenvalues of — A^i. 

1.5 Lipschitz nonlinearities 

Now, what happens when / is not smooth anymore, but, say, Lipschitz? In par¬ 
ticular, this is the scenario considered in the proof of the so called one dimen¬ 
sional Lazer-McKenna conjecture (i, [To]) by Costa, Figueiredo and Srikanth in 
[3]. We state the result, for the reader’s convenience. Let X = iL|,([0,7r]) be the 
Sobolev space of functions satisfying Dirichlet boundary conditions with square in- 
tegrable second derivatives. Recall that u —u" acting on X to F = L^([0, tt]) has 
eigenvalues Xk = k^, k = 1,2,... with corresponding eigenfunctions sin(fcx). Take 
/ : R —>■ R, a strictly convex smooth function / with asymptotic values a and b for 
its derivative /' satisfying 

a < 1, Xk = k"^ < b < {k + 1)"^ = A/c+i. 

Then the equation F{u) = —u" — f{u) = —tsinx, u{0) = u{tt) = 0, has exactly 2k 
solutions for < > 0 sufficiently large. 

The argument in [3] considers the nonlinearity / given by f'{x) = a or b, de¬ 
pending if cc < 0 or a: > 0. The related operator F(u) = —u" — f(u) now requires 
some care: it stops being differentiable everywhere and the usual differentiable nor¬ 
mal forms at regular points and folds break down. 

In this paper, we show that when / is merely Lipschitz, for appropriate condi¬ 
tions on the boundary of il, the operator F is still flat, in the sense that the global 
Lyapunov-Schmidt decomposition still holds, and sheets and fibers are still available 
as graphs of Lipschitz functions. A word of caution: piecewise linear nonlinearities 
may yield continua of points on which the map F takes a unique value. Such sets 
necessarily lie in a single fiber. More, the numerical analysis of the PDE F{u) = g 
presented in |4] is still valid, after minor modifications. 

We take this material to be an intermediate step towards a more geometric 
description of the operators of Hamilton-Jacobi-Bellman type, as studied by Felmer, 
Quaas and Sirakov in [7] . 

The authors gratefully acknowledge support from CAPES, CNPq and FAPERJ. 
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2 Some cautionary examples 

For a box ft C M", we consider again the differential equation 

F{u) = -Au-f{u) = g, M|an = 0 (2) 

with a Lipschitz nonlinearity f{x). Once / is piecewise linear, there may be whole 
(straight line) segments on the domain restricted to which F is actually constant. 
This happens already for the one-dimensional case. Let ifi be the (positive) ground 
state associated to eigenvalue Ai and take a < Xi < b. Suppose that Ran ipi = 
[0, M] and define f{x) to be continuous, with derivatives equal to a, Ai and b in the 
intervals [—oo, 0], [0, M] and [M, oo]: clearly, F{t(pi) = 0, for t G [0,1]. 

2.1 Nonlinearities / with derivatives taking two values, n = 1 

In a similar vein, we now provide examples of segments on which F is constant for 
the nonlinearity f{x) = ax or bx, for a; < 0 or x > 0 , with the property that, for 
special values of a and &, there are right hand sides g (5 = 0 is an example) with 
the property that F~^{g) contains a (straight) half-line of solutions. In particular, 
the map F : X = fl ilo(Al) —>■ is not proper and the equation has 

solutions which are not isolated. 



Figure 1: A half-line of solutions 


We begin with the one dimensional case, n = 1. Set If = / = [0,7r]. Split / 
into k closed intervals Ii,i = joined at their ends, of two different lengths, 

lodd i = P, leven i = «• In the figure, k = 4. The smallest eigenvalues for the 
operator u 1 —>■ —u" with Dirichlet conditions in an interval of sizes /3 and a. are 
respectively 

A^ = Q^ and A„ = (-)^ 
p a 

with positive (normalized) eigenfunctions tpp and ipa- We set b = \p and a = Xa and 
construct a solution ip by juxtaposing multiples of (pp and Pa as shown in the figure. 
On Ii , one may take ppp, for arbitrary p > 0. On I 2 , the (negative) multiple of ipo, is 
determined by matching the first derivative — recall that ip € X = U 

so Ip' is absolutely continuous. The procedure extends to the remaining intervals 
in a unique fashion. We have to make sure that the total length of the intervals Ii 
equals tt. Thus, for example, in the simplest case /c = 2, we must have 


(3 + a = t: 


I 

Vb 



For any value a G (1,4), there is a (unique) b, which turns out to be in (4, 00 ) which 
solves this equation. Said differently: any interval [a, b] containing A 2 = 4 for which 
a, b are not eigenvalues of the free problem admits a half-line of solutions of the 
equation above. For different numbers of intervals, one shows half-lines of solutions 
for any a G (Afc, Afc+i) and appropriate b G (Afc+i, 00 ). 
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Alas, the only situation for which this argument does not provide a half-line of 
solutions F“^(0) is a < Ai, the Ambrosetti-Prodi case. There are strong evidences 
that in this case there are no continua of solutions but we have no real 

proof. 

Clearly, one may replace 5 = 0 by any g defined piecewise on intervals A as 
functions in the range of rt 1 —>■ —u" — f{u) (Dirichlet conditions on li) acting on 
positive functions restricted to li. Thus, for k intervals, the set of such 5 is a 
vector subspace of L^{I) of codimension k. This construction ascertains that g is 
in the range of F (now considered in the full interval I), so that g = F{uo). By 
linearity on each interval li, adding an homogeneous solution ip S ^“^( 0 ) gives rise 
to uo + ')p € F~^{g). 

These ideas also suffice to prove that there is no nontrivial function in [0, tt] 
which is taken to 0 if o < Ai < 6 . 

2.2 The case of arbitrary n 

We now consider the case n = 2, the general situation being similar. We now have 
n = J X / (rectangles would work also): just separate variables and proceed. Let 
ipiy) as before, solving —ipyy — /('*/') = 0 ; where / is constructed from appropriate 
a and b, and let ip{x) > 0 be the ground state for Dirichlet conditions on I, so that 
= ‘p(x). The product ip{x,y) = ip{x)'ip{y) satisfies 

-Ax - Ay - fiA = -0 + (p{x){-Ay - /( 0 )) = 0 
so that 0 and its positive multiples solve 

- Uyy - f{u) = 0 , M|an = 0 , 


for f{x) = f{x) + X. 

3 Geometry of Lipschitz maps 

Set Y = L^(f2) with inner product {u,v)q = JqUv and norm ||u||o. Also let X = 
n iLg(D), with inner product {u,v} = (—Am , — Am)o and norm ||m||2. 

We always consider sets D for which —A : X C T —>■ F is a self-adjoint isomor¬ 
phism — we call such domains fl appropriate. Also the same operator should 
have as a core, i.e. it is essentially self-adjoint in this domain. From 

the spectral theorem, there is an orthonormal basis of (Dirichlet) eigenfunctions 
(fi € X, ||i^i||o = 1, satisfying —Aipi = Xi(pi. Eigenfunctions associated to different 
eigenvalues are orthogonal with respect to both inner products. Concretely, one 
might take D to be a convex set or require dfl to be C^’ ' ([n],is])- Notice that, 
from standard results in spectral theory, operators 

T-.XCY^Y, Tu = -Am - qu, 

for bounded real potentials q, are still self-adjoint with an orthonormal basis of 
eigenfunctions. 

We assume that the nonlinearity / : R —>■ K is Lipschitz, and /' takes values in 
an interval [a, b] with the property that the bounds a and b are not eigenvalues A^. 
For this part of the paper, we make no assumptions about convexity for /. Notice 
that a and b do not have to be the asymptotic values of /', a degree of freedom 
which is convenient for numerical analysis. 

For starters, F : X ^ Y given by F{u) = —Am — /(m) is a well defined map — 
it suffices to check that /(m) S Y. This follows from the easy lemma below. 
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Lemma 1 Say / : R —>■ R is Lipschitz. Then the map f : Y = LS‘{TI) —>■ Y given 
by f{u) = f o u is also well defined and Lipschitz with the same constant. 

Proof: Take first u, v continuous functions. Since / is M-Lipschitz, it is absolutely 
continuous, so that 

\f{u{x))\ = 1/(0) + u{x) [ f'{tu{x))dt I < |/(0)| + M\u{x)\,x G ff, 

Jo 

and, since fl is bounded, we have f{u) G Y. Similarly, applying the fundamental 
theorem of calculus to the function ip{t) = f{tu{x) + (1 — t)v{x)), one obtains 

\\f{u{x)) - fiv{x))\\o < [ \f{tu{x) + (1 - t)v{x))\dt ||u(x) - ?;(a;)||o. 

Jo 

Now take Cauchy sequences of continuous functions converging to arbitrary func¬ 
tions in Y: the estimates above extend to the required estimates. ■ 


3.1 The main result: Fy : Wx Wx is an isomorphism 

We now describe an orthogonal decomposition of X and Y. Take il C R” to be a 
bounded appropriate domain and let Ay = {Aijig/ be the set of eigenvalues Xi in 
(a, 6 ). The vertical subspaces Vx = Vy equal the invariant subspace associated to 
Ay and Vx C X,Vy C Y. The horizontal subspaces are Wx = C X and Wy = 
Vy C Y where orthogonality takes into account the (different) inner products in X 
andT. These induce orthogonal decompositions X = Wx®Vx, Y = Wy(BVy And 
corresponding orthogonal projections Py and Qy from Y to Wy and Vy. Finally, 
we consider affine horizontal subspaces in X, which are sets of the form x + Wx, 
for a fixed x, and affine vertical subspaces in Y, of the form y + Vy. 

We need a label for this construction: a nonlinearity / induces an /-decomposition 
X = Wx ®Vx, Y = Wy © Vy associated to bounds a and b. 

Theorem 1 Let LI be an appropriate domain, / : R —>■ R Lipschitz, Ran/' C [a, b], 
where a and b are not eigenvalues Xi. and X = Wx © Vx, Y = Wy © Vy be the 
I-decomposition specified above. For v G Vx, let : Wx —>■ Wy be the horizontal 
projection of the restriction of F to the affine subspace v-hWx, Fy(w) = PyF{w-\-v). 
Then Fy is a bi-Lipschitz homeomorphism. The Lipschitz constants for Fy and 
are independent ofv. 

Proof: For 7 = (a + b)/2, set f{x) = f{x) — jx. Then T : Wx —>■ Wy given by 
u —>■ —Am — 7 M is well defined and invertible, with eigenvalues Ay — 7 with j ^ /. 
Let Am — 7 be the eigenvalue of T of smallest absolute value: clearly, 

= |Am > |a- 7 | = 6 - 7 . 

For u = w -\- v,w € Wx, v G Vx , we have 

Fy{w) = Py[—A{w + m) — f{w + v)] = Tw — Pyf{w + v). 

The composition Fy o T~^ : Wy —>■ Wy is of the form / — Ky, where Ky(w) = 
Pyf{T~^w + v). We show that Ky : Wy —>■ Wy is a contraction with constant 
uniformly bounded away from 1. For w,w G Wy C L'^{Ll), 

IIA:^('u;) - Ar„(M;)||o < ||/(T“^m; + m) - f{T~^w + v)||o 

< (&-7)||T“^(m;-'u;)||o < ^ ||m;- w||o =c ||M;-ri;||o. 

I^m TI 
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Since b is not an eigenvalue Xj , the Lipschitz constant c is uniformly bounded away 
from 1. From the Banach contraction theorem, Fy : Wy Wy is a homeomorphism 
and standard estimates using the familiar formula 

(/ - Ky) ^ = / + Ky + Ky O Ky + Ky O Ky O Ky + . . . 

show that {Fy)~^ is Lipschitz, where the constant depends on c and not on v. I 

In particular, if [a, 6 ] does not contain eigenvalues Xi, the result above recovers 
the Dolph-Hammerstein theorem for Lipschitz nonlinearities (0, 0). 


3.2 Sheets and fibers make sense for Lipschitz nonlinearities 

We are ready to extend to the Lipschitz context the global Lyapunov-Schmidt de¬ 
composition which is known for the smooth case, from the works of Berger and 
Podolak and Smiley. Consider the following diagram. 

X = Wx © Vx -^ X = Wy © Vy 

•S>-^={Fy,Id) \ /' F=Fo'S>={Id,ij>) 

Y = Wy®Vy 


Thus the change of variables $ yields F{w,v) = F o ^{w,v) = {w,4>{w,v)), 
from which we will derive some convenient geometric properties. We first clarify a 
technicality: $ is indeed a global change of variables in the Lipschitz category. In 
Wx © Vx , we use the norm obtained by adding the norms in each coordinate. 



Figure 2: The change of variables 


Proposition 1 The map $ = {(Fy) ^,Id) : Y = Wy (BVy ^ X = Wx (BVx is a 
hi-Lipschitz homeomorphism. 

Proof: The invertibility of $ follows from the previous theorem. We use some 
elementary facts. The identity map Id : (VxjIMb) —>■ (Ly,||.||o) between normed 
spaces of the same hnite dimension is bi-Lipschitz. Also, since Wx and Vx are 
orthogonal (in and H'^), ||w ||2 + ||n ||2 < 2 ||w ± n ||2 for w G Wx and v G Vx- To 
show that is Lipschitz, take w + v,w + v G X. For appropriate constants C, C, 

+ ?;) - + f;)||o = \\Fy{w) - F{j(n})||o + ||u - v\\o 

< II - A{w -w)- Py{f{w + v) - f{w + I;))||o + C\\v - nib 
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< \\w - w \\2 + \\f{w + v)- f{w + i;)||o + C\\v - ■u ||2 < Cllw + d - (w + 0)||2, 


where the last inequality follows from Lemma 1. 

We obtain a Lipschitz estimate for $ = ((F„)“^, Id) : ILV © foy —>■ Wx © Vx- 
Take z + v,z + v€V = Wy © Vy. Then 

+ n) - $(z + 0)||2 < \\F-\v) - Fr\v)h + \\Fr\v) - F--\i;)\\, + ||^; - 

Again from finite dimensionality of Vx = Vy, there is an estimate of the form 
lit! — n ||2 < C\\v — i)||o. We also have \\F^^{v) — F^^{v )\\2 < C\\v — i;||o from the 
proof of the previous theorem. The first term is handled in a similar fashion. ■ 

The picture should help putting pieces together. Here, dimVA = dimVy = 1, 
as in the Ambrosetti-Prodi theorem: the convex span of Ran/' (/ is Lipschitz!) 
contains only the eigenvalue Ai. The map F takes an affine horizontal subspace 
V + Wx to a sheet, and the inverse of the vertical affine subspace 5 + Vy is a fiber, 
which crosses Wx at u;( 0 ) and v + Wx at w + w{v), in the notation of the proof 
above. Clearly, sheet and fiber are graphs, as stated above. The change of variables 
$ preserves horizontal affine subspaces and F preserves affine vertical subspaces. 

We are ready to prove the fundamental geometric property of such F : X ^ Y: 
there are uniformly flat sheets and uniformly steep fibers. 


Proposition 2 Let F : X = Wx(BVx Y = Wy©Vy with the hypothesis given in 
the beginning of the section. The image of each horizontal affine space v + Wx C X 
under F is the graph of a Lipschitz function : Wy -P- Vy . Similarly, the inverse 
of each vertical affine subspace g + Vy C Y under F is the graph of a Lipschitz 
function ag : Vx —>■ Wx- The Lipschitz constant can be taken to he the same, for 
all V G Vx, g € fVy. 

Proof: We prove the result for fibers F~^{g + Vy): the statement for sheets 

F{v + Wx) is easier. Clearly, F~^{g + Vy) C g + Vy, which is taken by the 
change of variables $ to a set of the form ag = {{Fy)~^{g) + v,v £ Vx} C X. 
From the theorem, for every v G Vx, there is a unique w(y) G Wx for which 
PyF{w{v) + v) =g — thus F“^((/ + Vy) = {iFy)~^(g)+v,v £ Vx}- Said differently, 
w{v) -£ V £ ag'. the set {(w{v), v),v £ Vx} C A is the graph of a Lipschitz map. 

The uniformity (on g) of the Lipschitz constant of the maps v 1 —>■ w{v) is respon¬ 
sible for the uniform steepness of the fibers. ■ 

In opposition to the arguments in 0 , m and [3] for the smooth case, the 
geometric statements follow without recourse to implicit function theorems. Notice 
also that the uniform flatness of sheets and steepness of fibers are a counterpart to 
(differential) transversality between fibers and horizontal affine spaces in the domain 
and between sheets and vertical affine spaces in the counterdomain. 

The restriction of F to horizontal affine subspaces is injective but the restriction 
to fibers a is not. In particular, in the standard Ambrosetti-Prodi case, F restricted 
to each fiber is simply the map x G R 1 —>■ —x^ £ M, after global changes of variables. 
The theorem becomes evident from this fact, first proved in [2]. 

Vertical lines may be taken by F to the horizontal plane, indicating yet another 
relevant transversality property of the fibers. To see this, take H = [—7r/2,7r/2] and 
F{u) = —u" — f{u), so that A 2 = 4 and the corresponding eigenvector (p 2 is odd. 
Set a = 3 and b = 5, split X = Wx © {^ 2 ) and take f{x) = e{x) + 4x, where e{x) 
is even (convexity is not necessary!) and Ran/' = (a, b). By symmetry, we have 


= [ {-tA(p2 

JQ, 


e{t(p2) — 4t 1 P 2 ) T2 = 0. 


( F{tip2) , T2 ) 


4 Geometry and numerics 

The statements in the previous section are exactly what we need to mimic the 
numerical algorithms in [3] for solving F{u) = g, i.e., the differential equation (2). 
This section emphasizes the points where some alterations are needed, but most 
details common to the smooth and Lipschitz scenarios, which are provided in [4], 
will not be presented. 

4.1 Finding the right fiber 

To solve F{u) = g, first find any point u = w+v € ag, the fiber associated to g. Said 
differently, find u so that PyF{u) = Pyg- In order to do this, notice that, from the 
results of the previous section, each fiber ag intersects each horizontal affine space 
V + Wx at a single point. Thus each horizontal point Pyg G Wy corresponds to a 
unique w{g) G Wx for which Fy(w{g)) = Pyg. 

Said differently, each bi-Lipschitz map : Wx — Wy takes a point 'w{g) in the 
fiber ag to a point Pyg, which is the only point in the vertical affine space g + Vy 
in the horizontal subspace Wy. In words: there is a bi-Lipschitz map between the 
set of all fibers (represented by points in Wx) to the set of vertical affine spaces 
(represented by points in Wy). 

Now the good (numerical) news: to invert each map Fy, simply invert the (bi- 
Lipschitz) homeomorphism Fy o = I — K, where K is a contraction, as shown 
in the proof of Theorem I! So the approximation of F~^{PyG) is amenable to 
standard numerical algorithms. In the differentiable case, we could do somehow 
better: we could invert by continuation where local steps are given by Newton 
iterations. In the strictly Lipschitz context, we lose quadratic convergence — some 
acceleration techniques are still available, but we provide no details. 

4.2 Moving along a fiber 

Once a point in ag is identified, we need to learn to walk along the fiber, or more 
precisely, we have to compute the point in the fiber with a given height v G Vx • This 
is turn is the main piece of information for a finite dimensional inversion algorithm 
for the restriction F : ag ^ g + Vy. This is done by continuation starting from a 
given point in the fiber u G ag. 

Clearly small perturbations of a point ui G ag leave the fiber. But the previous 
algorithm — more precisely, the inversion of each map Fy — makes it possible to 
change ui = wi + vi to a point u = w + V 2 and then using ft as a starting point to 
solve Fy.^(w 2 ) = Pyg, giving rise to a point U2 = W2 + V2 G ag. 

4.3 Stability of the decomposition, a numerical necessity 

The uniformity on the flatness of sheets and the steepness of fibers has a relevant 
consequence for numerics, which has been detailed in [4]. In a nutshell, whatever 
algorithm we use requires an approximation for the projection Py : Y ^ Y. This 
is usually accomplished by computing (approximate) eigenfunctions associated to 
the eigenvalues of —Ad in the interval {a,b). A concrete possibility is to approx¬ 
imate functions by finite elements. The upshot is that the numerics handles an 
approximation Vx = Vy of Vx = Vy and its orthogonal complement. However, 
from the uniformities, such spaces, for sufficiently close approximations, still in¬ 
duce global Lyapunov-Schmidt decompositions giving rise to sheets and fibers, for 
which the algorithms described in the previous paragraphs hold and provide robust 
approximations to the real answers. 
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Figure 3: Irrelevant perturbations for numerical purposes 


In the figure, we sketch fibers and their perturbations; there is an analogous 
reasoning for sheets. In the figure, Vy changes slightly the affine vertical subspace 
through the point g, which in turn, when inverted, gives rise to a slightly different 
fiber dg. Still, dg is a Lipschitz graph of a function from Vx to Wx = (Vx)^- The 
(geometric) thing to notice is the fact that it is the steepness of the fibers which 
ascertain the robustness of the global Lyapunov-Schmidt decomposition. 
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